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ABSTRACT 

Wc characterise the stellar masses and star formation rates in a sample of -^40,000 
spectroscopically confirmed UV luminous galaxies at 0.3<z<1.0 selected from within 
the WiggleZ Dark Energy Survey. In particular, we match this UV bright population to 
wide-field infrared surveys such as the near infrared (NIR) UKIDSS Large Area Survey 
(LAS) and the mid infrared WISE All-Sky Survey. We find that ~30 per cent of the 
UV-luminous WiggleZ galaxies, corresponding to the brightest and reddest subset, are 
detected at >5cr in the UKIDSS-LAS at all redshifts. An even more luminous subset of 
15 per cent are also detected in the WISE 3.4 and 4.6/im bands. In addition, 22 of the 
WiggleZ galaxies are extremely luminous at 12 and 22/xm and have colours consistent 
with being star formation dominated. We compute stellar masses for this very large 
sample of extremely blue galaxies and quantify the sensitivity of the stellar mass esti- 
mates to various assumptions made during the SED fitting. The median stellar masses 
are logio(M*/MQ)=9.6±0.7, 10.2±0.5 and 10.4±0.4 for the IR-undetected, UKIDSS 
detected and \JKIDSS+ WISE detected galaxies respectively. We demonstrate that 
the inclusion of NIR photometry can lead to tighter constraints on the stellar masses 
by bringing down the upper bound on the stellar mass estimate. The mass estimates 
are found to be most sensitive to the inclusion of secondary bursts of star formation 
as well as changes in the stellar population synthesis models, both of which can lead 
to median discrepancies of the order of 0.3 dex in the stellar masses. We conclude that 
even for these extremely blue galaxies, different SED fitting codes therefore produce 
extremely robust stellar mass estimates. We find however, that the best-fit M/Lk 
is significantly lower than that predicted by simple optical colour-based estimators 
for many of the WiggleZ galaxies. The simple colour-based estimator over-predicts 
M/Lk by ~0.4 dex on average. The effect is more pronounced for bluer galaxies with 
younger best-fit ages. The WiggleZ galaxies have star formation rates of 3-10 M0yr~^ 
and mostly lie at the upper end of the main sequence of star-forming galaxies at these 
redshifts. Their rest-frame UV luminosities and stellar masses are comparable to both 
local compact UV-luminous galaxies as well as Lyman break galaxies at z 2 — 3. The 
stellar masses from this paper will be made publicly available with the next WiggleZ 
data release. 

Key words: galaxies: stellar content, galaxies: formation, galaxies: evolution 



1 INTRODUCTION 

In recent years, large area spectroscopic surveys of both pas- 
sive galaxies such as the Luminous Red Galaxy (LRG) pop- 
ulation, as well as emission line galaxies (ELGs), have been 
successfully used to place accurate constra i nts on cosmo- 
logica l models (e.g [Blake et al.l |2007|, 120081 : iPercival et"al] 
I2OIOI : iBlake et al.] I2OIII : iReid et all I2OI2I '). At the same 
time, deeper, smaller area spectroscopic surveys such as 
DEEP2 and VVDS as well as multi- wavelength photomet- 
ric surveys such as C OMBO-17 (iBell et al.ll2004ll. CFHTL S 
lArnouts et al.1 1200/1 ) and COSMOS jllbert et al.1 l2009t ). 
have enabled detailed studies of galaxy formation and evo- 
lution and constraints on the global properties such as the 
spectral energy distributions (SEDs), ages, stellar masses 
and star formation histories (SFHs) of galaxies out to z~2. 
The advent of very large area photometric surveys at multi- 
ple wavelengths now offers us the possibility of constraining 
the global properties of the large spectroscopic samples that 
have been assembled for cosmology. Although the large area 
surveys are by their nature very shallow, they contain huge 
numbers of galaxies at z<l, therefore enabling a statistically 
robust census of galaxy properties at these redshifts. 

Star-forming galaxies at the main epoch of galaxy for- 
mation at z~l-3 have been selected in many different 



ways. The most luminous starbursts at these epochs are of- 
ten selected at long w avelengths such as the far infrared 
and submillimeter (e.g . Small et al.ll2002l: llvison et al.ll2002l: 
Chapman et al.1 12003 : iMagnelh et al.1 I2OI0I : iBanerii et"aLl 



-ap: 

201lf ) while more modest star-forming galaxies have been 



targeted using optical colour - cuts s uch as the Lyman break 
selection (e.g. IPettini et al.l l200lh. the B M/BX method 
(|Adelberger et al.l |2004; ISteidel et al] \2004 ) and the BzK 
technique (|Daddi et al.1 120041 '). The availability of UV data 
from the Galaxy Evolution Explorer (GALEX) has also al- 
lowed comprehensive s tudies of UV-luminou s galaxies both 
in the loc al Universe |Heckman et al.ll2005|) and at hi gher 
redshifts (|Burgarella et al.ll2006l : iHaberzettl et al.ll2011^ . 

In this work, we study the physical properties of a pop- 
ulation of UV-luminous emission line gala xies selected from 
withi n the WiggleZ Dark Energy Survey (jPrinkwater et all 
l2010t ). Although primarily designed as a cosmology survey 
targeting blue emission line galaxies at intermediate red- 
shifts of z~0.7, the dataset contains ~215,000 spectroscop- 
ically confirmed highly star-forming galaxies that form a 
very large statistical sample that can also be exp loited for 
galaxy evolution studies (e.g. lWisnioski et al.ll201ll . Jurek et 
al. in preparation). In particular, the redshift range of the 
WiggleZ survey coupled with the UV selection means this 
sample is particularly useful for bridging the gap between 
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analogously selected local galaxies and UV-luminous galax- 
ies at z>l. While the large sample size makes targeted multi- 
wavelength follow-up of the majority of WiggleZ galaxies un- 
feasible, one can use existing multi-wavelength photometric 
datasets to better characterise this population. 

Identifying synergies between existing multi-wavelength 
photometric surveys and large redshift surveys like WiggleZ, 
is important for several reasons. Firstly, spectroscopic sam- 
ples such as the WiggleZ sample serve as important calibra- 
tion sets for large photometric surveys. As an example, the 
ongoing Dark Energy Survey (DES) will detect 300 million 
galaxies out to z~2 and will overlap the all-sky near infrared 
(NIR) VISTA Hemisphere Survey (VHS). The combination 
of optical and NIR data will allow photometric redshifts 
and SED fit parameters to be constrained for a v ery large 
number of galaxies out to z~2 jBanerii et al.|[2008l ). Colour 
selected spectroscopic samples like WiggleZ may be used 
Eis training sets but will by their nature be incomplete in 
certain regions of parameter space sampled by flux-limited 
photometric surveys such as DES and VHS. Quantifying this 
incompleteness in terms of the physical properties of the 
galaxies - i.e. understanding the types of galaxies that con- 
stitute currently available spectroscopic samples, is a useful 
exercise in order to better calibrate redshifts from photomet- 
ric sur veys. Next generatio n spectrosc opic surveys like Big 
BOSS (ISchlegel et al.ll201l|V 4M0 ST (|de Jong et al.l 120121 ) 
and DESpec (|Abdalla et al.ll2012l ). will in turn use the pho- 
tometric surveys as the basis for target selection. Once again, 
understanding the physical properties of existing spectro- 
scopically confirmed galaxies within these surveys, will help 
to design colour selection algorithms for new populations in 
the future. 

In this work, we characterise the stellar masses and 
star formation histories of a large sample of ~40,000 galax- 
ies at 0.3< z <1.0, corresponding to a single field in the 
WiggleZ survey, and focus in particular on the subset of 
^12,000 of these galaxies that are also matched to the NIR 
UKIDSS Large Area Survey (ULAS). The advantage of long- 
wavelength data in the infrared is that it is less sensitive to 
the dust extinction in galaxies than the UV and optical. 
This, combined with the fact that older, more massive stars 
are brighter in the infrared, means that the infrared data is 
expected to give a more unbiased estimate of the total stel- 
lar mass in galaxi es than the UV/optical (|Bell et al.ll2003l : 
iDrorv et al. I l2004h . Recent studies have however claimed 
that the NIR actually leads to poorer stellar mass esti- 
mates for a sample of intermediate redshift gala xies from the 
Gala xy and Mass Assembly (GAMA) survey (|Tavlor et al.l 
I2OIII ). This worsening in the stellar mass estimates is at- 
tributed to uncertainties in the stellar population synthesis 
models at these wavelengths and/or uncertainties in cali- 
brating the photometry between the optical and NIR sur- 
veys. The role of NIR photometry in constraining the stellar 
masses of galaxies therefore still remains open to debate. 

The stellar mass estimates themselves are useful for cos- 
mology given that the clustering strengths of galaxy sam- 
ples split by ste llar mass are expected to be different - (e.g. 
ICoil et al.l 1200^ 1. Our aim in the current work is therefore 
also to test the robustness of stellar mass estimates from 
SED fitting codes, to various assumptions made during the 
fitting process. The WiggleZ galaxies are sele cted to be ex- 
tremely blue galaxies (|Drinkwater et a"l]|2010l ). and as such. 



represent a sample where SED fitting is likely to be the most 
problematic. By testing SED fitting codes on this sample, 
we can be reasonably confident that the codes can be ap- 
plied to older, redder and more massive galaxies, which are 
likely to have less complex star formation histories. 

Throughout this paper we assume a fiat ACDM cos- 
mology with h—0.7. All magnitudes are on the AB sys- 
tem where the UKIDSS photometry has b een converted 
to the AB system using the conversions in iHewett et al.l 
(|2006h : Y=-H0.633, J=-h0.937, H=-(-1.376 and K=-Hl.897. 
The Vega to AB conversions in the WISE bands are as- 
sumed to be Wl=+2.683, W 2=-H3.319, W3=-F5.242 and 
W4=+6.404 (|Cutri et al.ll2012l ). 



2 DATA 

We begin by describing the spectroscopic and photometric 
catalogues that are used in this work. 



2.1 Wiggle-Z Dark Energy Survey 

Spectroscopic data for the UV luminous emission line 
galaxies is taken from the WiggleZ Dark Energy Survey 
(fPrinkwater et al.|[201oh . The survey has assembled reliable 
redshifts for 219,682 galaxies in seven different fields. The 
spectroscopic targets a re selected using ult raviolet data from 
the GALEX satellite (|Martin et al.1 120051 ') and optical data 
from the Sloan Digital Sky Su rvey (SDSS) in the north 
(|Adelman-McCarthv et al] |2006| ) and the Canada France 
Hawaii Telescope Red Sequence Cluster Survey in the south 
(|Yee et al.|[2007l ). A series of magnitude and colour cuts 
are applied to the data to preferentially select star-forming 
galaxies with bright emission lines which are then targetted 
using the AAOmega spectrog raph on the Anglo- Australian 
Telescope (|Sharp et al.l l2Q0q ). Full detai ls of the spectro- 
scopic target selection can be found in iDrinkwater et al] 
(|2010. ) but results in a complicated selection function over 
th e total survey area . This selection function is computed 
m iBlake et al.l (|2010D where the different sources of incom- 
pleteness are fully detailed. We note in particular that the 
complicated colour selection of WiggleZ targets means the 
sample is incomplete in some regions of redshift, stellar mass 
and star formation rate parameter space. 

In this paper, we work with the WiggleZ data in the 
15hr field only with 209<RA<231 and -3.2<DEC<7.2 in 
order to ensure overlap with currently available infrared 
datasets from UKIDSS and WISE. This field contains 46,144 
galaxies in total down to a fiux limit of NUV<22.8. All pho- 
tometric catalogues from the WiggleZ survey contain the de- 
reddened galaxy magnitudes and this sample only contains 
sources with redshift quality between 3 and 5 which corre- 
sponds to reliable redshift estimates. We select only galaxies 
at 0.3<z<1.0 for this study which constitutes the bulk of the 
WiggleZ population. At higher redshifts, the quality of the 
redshift estimate becomes increasingly unreliable and many 
of the z > 1 sources with reliable redshift measurements are 
AGN. Although these AGN may contain emission related 
to star formation, the SED fitting codes used in this work 
do not allow us to disentangle the contributions of the two. 
We are primarily interested in the star-formation dominated 
WiggleZ galaxies in this work and so we restrict our redshift 



© 0000 RAS, MNRAS 000, 000-000 



4 



0.3 < z < 1 .0 
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Figure 1. Observed r-band magnitude versus (g — r) colour for 
all WiggleZ galaxies in the 15hr field at 0.3< z <1.0, compared 
to star-forming galaxies from the DEEP2 survey and Luminous 
Red Galaxies from the BOSS survey, both at similar redshifts to 
WiggleZ. The contours represent the density of points computed 
using a kernel density estimator over a grid traversing the pa- 
rameter space covered by the samples. The WiggleZ galaxies are 
brighter than the DEEP2 sample but slightly bluer than these 
DEEP2 galaxies which represent normal star-forming galaxies at 
these redshifts. 



range to a regime where the AGN make up an insignifi- 
cant proportion of the population. Our final sample of Wig- 
gleZ galaxies therefore totals 39,701 sources at 0.3<z<1.0. 
The WiggleZ magnitudes used throughout this paper corre- 
spond to the de-reddened model magnitudes from SDSS and 
GALEX. The WiggleZ galaxies are typically ~5(7 detections 
in the GALEX NUV band and >10a detections in the SDSS 
r-band. 

In order to better illustrate the types of galaxies se- 
lected using the WiggleZ colour cuts, in Figure [1] we show 
the observed colour-magnitude distribution of the WiggleZ 
gala^xies. This is compared to ^19,000 galaxies from the 
DEEP2 survey Data Release 4 (| Newman et aLlbOld ) over 
the same redshift range. These DEEP2 galaxies represent 
typical star- forming galaxies at these redshifts. We also com- 
pare the distribution to ~180,000 Luminous Red Galaxies 
(LRGs) over the same redshift range selecte d from within 
the SDSS BOSS survey (jMaraston et al.ll2012l '). As expected, 
the WiggleZ galaxies are considerably bluer and fainter than 
the LRGs that make up the red sequence at these redshifts. 
However, despite the brighter flux limit of WiggleZ com- 
pared to DEEP2, at a fixed r-band magnitude, the WiggleZ 
selection targets sources that are also bluer than the typ- 
ical blue cloud galaxies that make up the DEEP2 sample. 
In other words, the WiggleZ selection is isolating the most 
extreme end of the blue galaxy population at any given lu- 
minosity. 

Having described the properties of the WiggleZ sam- 
ple, we now move on to considering synergies between this 
sample and wide-field infrared (IR) imaging surveys. 



2.2 UKIDSS Large Area Survey 

The UKIDSS Large Area Survey (|Lawrence et al.l 120071 ') is 
the current largest NIR survey and has obtained imaging 
over -^3200 deg^ of the northern sky in the Y,J,H and 
K bands. The survey is being carried out using the Wide 
Field Camera (WFCAM) on the 3.8m UK infrared Tele- 
scope (UKIRT). UKIDSS, which began in 2005, is the suc- 
cessor to the Two Micron All Sky Survey (2MASS) and is the 
NIR counterpart to the Sloan Digital Sky Survey (SDSS). 
We use Data Release 9 of the UKIDSS LAS (ULASDR9) 
in this work, which reaches nominal 5a depths of F=20.8, 
J=20.5, H=20.2 and K=20.l. Throughout this work, we 
use the Petrosian magnitudes in the UKIDSS catalogues as 
the UKIDSS catalogues do not include model magnitudes. 
These Petrosian magnitudes serve as a reasonable estimate 
of the total NIR flux of the galaxy. Differences between the 
SDSS Petrosian magnitudes and SDSS model magnitudes o f 
galaxies are of the order of 15 per cent (|Banerii et al.|[201ol V 
We match the WiggleZ galaxies in the 15hr fleld to 
UKIDSS using a matching radius of 2". The median sep- 
aration between the WiggleZ and UKIDSS sources is ~0.3". 
11,919 of the original 39,701 galaxies are detected at >5a in 
at least one of the UKIDSS bands corresponding to 30 per 
cent of the WiggleZ sample. The redshift distribution for this 
UKIDSS detected sub-sample is very similar to that of the 
entire population at 0.3< z <1.0 with a peak at Z'^0.7. The 
UKIDSS sub-sample is compared to those WiggleZ galaxies 
not matched to a NIR source in Figure [S] We flnd, as ex- 
pected that the fraction of very blue galaxies detected in the 
NIR is very low and increases as we go to redder colours. 
The UKIDSS detected galaxies are also brighter than those 
that are undetected. The median r-band magnitude of the 
UKIDSS detected galaxies is 21.0 versus 21.8 for the galax- 
ies undetected in UKIDSS. The median [g — r) colour is also 
0.2 magnitudes redder for the UKIDSS detected galaxies. 



2.3 Wide-Field Infrared Survey Explorer 

The Wide-Field Infrared Survey Explorer {WISE; 
IWright et all |2010| ) has conducted an all-sky survey in 
four passbands - 3.4, 4.6, 12 and 22fim with 5a depths 
of >19.1, >18.8, >16.4 and >14.5. By its nature, WISE 
contains many different populations of astrophysical sources 
including planetary debris disks, populations of cool low- 
mass stars and ultraluminous infrared galaxies and active 
galactic nuclei. Although it is very shallow at 12 and 22/im, 
the 3.4 and 4.6/im depths are reasonably well-matched to 
the UKIDSS-LAS so it is interesting to ask what fraction of 
very blue star-forming galaxies such as the WiggleZ sample, 
are detected at these wavelengths in a relatively shallow 
all-sky IR survey. 

We match our sample of 39,701 WiggleZ galaxies at 
0.3<z<1.0 to the WISE All-Sky Release using a matching 
radius of 4". The angular resolution at 3.4 and 4.6/.im is 
~6". We select only those galaxies that are detected at 
S/N>5 at 3.4/im and S/N>3 at 4.6/im and have WISE 
colours of [3.4^,„— 4.6^m] <0.8, consistent with these galax- 
ies n ot having a significant AGN component (|Assef et al.l 
I2OIOI ). At the WISE wavelengths, the presence of an AGN 
can significantly affect the galaxy colours and the SED fit- 
ting codes used in this work do not allow us to fit for this 
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AGN component. The median separation between the WISE 
and WiggleZ sources is ~0.5" and as expected, larger than 
that between the UKIDSS and WiggleZ sources. Note that 
using a smaller matching radius of 2" rather than 4", de- 
creases the number of galaxies by only ~10% and does not 
affect any of our results. We use the W ISE magnitudes ob - 
tained from profile fitting as detailed in lCutri et al.l (|2012h . 
which mitigates the effects of source confusion and source 
blending which can be significant at the typical fluxes of 
our sample. At the redshifts of the WiggleZ sample, we 
note that most of the galaxies appear unresolved in the 
WISE images. More than 90% of the WISE matches are 
also in UKIDSS. Visual inspection of those that are not 
shows them to mostly be blended sources where two nearby 
sources that are resolved separately in the UKIDSS im- 
ages, are blended together in WISE. The remaining WISE 
only identiflcations correspond to small patches of sky with 
no UKIDSS coverage. We therefore only restrict our WISE 
sample to those sources that are also present in UKIDSS. 
The VF/Si^-I-UKIDSS sub-sample contains 6117 galaxies at 
0.3< z < 1.0. Their colour-magnitude distribution can also 
be seen in Figure [2] The WISE sample is brighter and 
slightly redder than the UKIDSS sample and corresponds 
to the brightest and reddest end of the WiggleZ popula- 
tion. The median r-band magnitude for the WISE detected 
galaxies is 20.7 and these gala^xies are 0.05 magnitudes red- 
der on average than the UKIDSS detected galaxies. Once 
again we find very little dependence of the matched fraction 
on the source redshift with the WISE sample also peaking 
at z~0.7. 

We also note that, despite the very shallow fiux limit in 
the WISE 12 and 22/im band, there are 78 WiggleZ galax- 
ies that are detected at >5a at these longer wavelengths. 
Thirty-seven of these also have [3.4^^— 4.6^m] <0.8, consis- 
tent with these galaxies not having a significant AGN com- 
ponent. A further sub-sample of 24 of these sati sfy the more 
conse rvative colour cut of [3.4^™— 4.6^m] <0.7 jStern et al.l 
|2012| ) and can reasonably be assumed to be star-formation 
dominated. We characterise the SEDs of these UV-luminous 
mid-IR bright sources in Section [S] The remaining Wig- 
gleZ galaxies with mid infrared emission can reasonably be 
thought to be AGN dominated and are not considered fur- 
ther in this work. 



22 22.5 




Figure 2. Obscrvcd-framc (g — r) colour magnitude diagram of 
infrared detected and infrared undetected WiggleZ galaxies show- 
ing that the UKIDSS-LAS and WISE detected sub-samples con- 
stitute the brighter, redder end of the WiggleZ population. The 
greyscale represents the density of points. The redshift distribu- 
tions for all three samples are very similar and peak at z = 0.7. 



3 SPECTRAL ENERGY DISTRIBUTION 
(SED) FITTING 

A key aim of this work is to derive SED fit parameters and in 
particular stellar masses, for the various sub-samples of UV- 
luminous galaxies from WiggleZ and to test the robustness 
of these parameters to assumptions made during the fitting 
process. In particular, we would like to test how these various 
assumptions affect the stellar mass estimates as a function 
of the UV-luminosity. We therefore begin by describing the 
various SED fitting codes that are used in this work and 
highlight similarities and differences between them. 



3.1 Fitting and Analysis of Spectral Templates 
(FAST) 

The FAST code is fully described in lKriek et all (|2009l ) and 
is designed to fit a range of different stellar population mod- 
els to fluxes of galaxies in the optical and NIR bands. It has 
been succesfuUy applied to infer the global properties of NIR 
selected samples of gal axies out to z^2 in the NEWFIRM 
medium band surve v (IKriek et al. I2OIOI). Cu rrently FAST 
allows fittii igof bothlBruzual fc CharlotI (| 2003h - BC03 here- 
after - and liviaraston ( 2005h models at a range of metallici- 
ties and using different initial mass functions (IMFs) and 
single component star formation histories (SFH) such as 
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exponentially decaying SFHs, delayed exponentially decay- 
ing SFHs and truncated SFH s. The BC03 m ode ls can be 
used in conjunction with the ISalpeteJ (|l955h or IChabrieil 
l|2003l ) IMF, while the Maraston mod el lib raries are con - 
structed using both the ISalpeted (|l955h and lKroupal (|200ll ) 
IMFs. Interste l lar du st extinction is accounted for using the 
ICalzetti et al] l|200Ct ) extinction law and is allowed to vary 
over a range of Ay- 

The advantage of FAST is that once a grid of model 
fluxes has been computed, the SED fitting is relatively quick 
as the code uses a Monte Carlo sampling of the model pa- 
rameter space to determine the best-fit parameters for each 
observed galaxy. This Monte Carlo sampling also allows 
FAST to calibrate the confidence intervals for each SED 
fit parameter that is estimated. Errors on the SED fit pa- 
rameters are calibrated using fOO Monte Carlo simulations. 
Therefore, the code provides realistic errors on the SED fit 
parameters that take into account both the uncertainties in 
the broadband fluxes used for the fitting as well as the un- 
certaintie s in the SED models v ia the use of a template error 
function (|Brammer et al.ll2008l ). These errors along with the 
reduced metric output by the code, can be used to assess 
the quality of the fits for different choices of input parame- 
ters. 

FAST is readily applicable to the large sample sizes cur- 
rently being assembled in cosmological volume surveys due 
to the speed with which it computes the best-fit SEDs. It 
does not however currently allow the inclusion of episodic 
bursts of star formation and the galaxy templates included 
do not contain emission lines which may be important for 
our sample of UV luminous gala:xies. FAST also does not 
include inverse-r or exponentially increasing star-formation 
histories which m ay be more representat ive of actively star- 
forming galaxies ( Maraston et al.l I2OIC1I ). Finally, we note 
that the SED models used by FAST do not account for 
dust emission in the rest-frame infrared and assume that the 
SEDs of the galaxies over the rest-frame wavelength range 
probed by our data, are dominated by starlight. Recent stud- 
ies have however noted the presence of excess emission at 
near infrared wavelengths that cannot be accounted for sim- 
ply by the stellar continuum. This excess emission is best 
modelled as an additional greybody compo nent with a tem- 
perat ure of between 750K and f 200K (e.g. Ida Cunha et al.l 
|2008| ). Some studies have found that the emission corre- 
lates with star formation and have therefore attributed it 
to emission from circumstel lar disks around massive stars 
(|Mentuch et all |2009| . I2OI0I '). However, as noted in these 
studies, this excess dust emission above the stellar con- 
tinuum only starts to contribute at rest-frame wavelengths 
above ~2/^m. 



3.2 MAGPHYS 



The MAGPHYS code (|da Cunha et al.1 12OO8I . l20ffl ) is de- 
signed to consistently treat the combined UV, optical and 
infrared emission from galaxies. In this code, any attenu- 
ation of starlight by dust at bluer wavelengths appears as 
reprocessed thermal emission in the infrared where the dust 
is treated in a physically consistent way bet ween the differ- 
ent wavelengths. The code uses the updated iBruzua] (|2007l ) 
model SEDs which include a new prescription for the TP- 
AGB stars to better reproduce the NIR colours of interme- 



diate age stellar populations. An advantage over FAST is 
that the code also allows for random bursts of star forma- 
tion to be added to the simple exponentially decaying SFH 
with bursts occurring with equal probability at all times 
since the formation redshift. The burst probability is set so 
that up to 50% of the galaxies in the model library have 
experienced a burst within the last 2 Gyr. The fraction of 
stellar mass formed in the burst ranges between 0.03 and 
4x that formed through continuous star formation with a 
characteristic timescale r. A consequence of the facility to 
incorporate more complex star formation histories as well 
as the facility to compute the IR emission in the galaxies in 
a self consistent way, is that the code is also slower. MAG- 
PHYS compares the observed photometry of galaxies to a 
total of ~66f million models and the priors imposed on the 
parameters are deliberately not overly restrictive so as to en- 
sure that the entire multi-dimensional observational space is 
reasonably well sampled. However this very large parameter 
space also means that the code is not as easily applicable to 
very large samples of galaxies such as the one in this paper. 

As the main advantage of this code over the others, is 
in the consistent modelling of the UV to infrared emission in 
galaxies, we only use MAGPHYS to characterise the prop- 
erties of the small subset of WiggleZ galaxies that are bright 
at 12 and 22fiia as well as to carry out independent checks 
on some of our results obtained using the other codes, using 
smaller random sub-samples of the WiggleZ galaxies. 



3.3 KG04 Code 

We also use a non-public code developed by one of us (KG) 
in the Perl Data Languagfl which has been used in se veral 
papers (e.g. IClazebrook et al.|[2004l : iBaldrv et al.ll2008l '). We 
will refer to this code as th e 'KG04 code'. This code uses the 
PEG ASE.2 SPS models (|Fioc fc Rocca-Volmerangj Il997l . 
If 999t ) and incorporates two-component star-formation histo- 
ries with a primary (long term) SFH plus a burst component. 
The primary component is represented by an exponentially 
decaying burst (with values r = 0.1,0.2,0.5,1,2,4,8,500 
Gyrjfl and the burst component is represented as a rapid 
exponential SFH of r = 100 Myr. The burst is allowed to 
occur over the full range of a galaxy SFH and is allowed 
to contribute a range of final stellar mass ratios from zero 
up to twice the long-SFH mass (evaluated at 13 Gyr). One 
important difference from the MAGPHYS code (see Section 
I3.2|l is that the burst fraction and times are not random but 
distributed across a full grid range. Dust is allowed to vary 
over a range of Ay with a Calzetti law, metallicity is held 
fixed with time in the PEGASE.2 models but is allowed to 
vary in the fitting. Nebular emission lines can optionally be 
included altho ugh we note that the line ratios a re held fixed 
in the models (|Fioc fc Rocca- Volmerangelll999l ). 

In a similar fashion to FAST the computation is sped 
up by pre-computing grids of model photometry values at 
each redshift of interest, and mass values and errors of a par- 
ticular galaxy are then determined by a fast minumum-x^ 
lookup of the photometry with a Monte-Carlo realisations 
of the errors. The final mass and error values are the mean 



^ http://pdl.perl.org 

r = SOOGyr is intended to approximate constant SFR models. 
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and standard deviation of five Monte-Carlo realisations. The 
grids in the various parameters are stepped in a pseudo- 
logarithmic fashion because of (i) the optimisation of the 
sampling of each parameter with respect to physical varia- 
tions and (ii) the constrained values of certain parameters 
such as the fixed metallicity values of the PEGASE.2 mod- 
els. A total of 2,427,480 possible model SEDs are considered 
for each galaxy in the fitting (with some being automati- 
cally excluded by the age of the Universe at a given red- 
shift constraint). Normally t he code has been used with the 
iBaldrv fc Glazebrookl (|2003l ) IMF (BG03 hereafter), how- 
ever it can be adapted to any IMF as is done in this paper. 

Stellar masses are the masses locked up in luminous 
stars at the best fit time but exclude mass contributions 
from non-luminous stellar remnants: white dwarfs, neutron 
stars and black holes. 



4 RESULT OF SED FITTING 

In this section, we test the robustness of stellar mass esti- 
mates for the WiggleZ galaxies to various assumptions made 
during the SED fitting procedure and consider in particu- 
lar the dependence of the mass estimates on the rest-frame 
FUV luminosity. 

4.1 The Effect of Photometric Bands 

We begin by assessing the role of infrared photometry in con- 
straining the stellar masses of the WiggleZ galaxies. There 
are several reasons why infrared data is expected to help 
in constraining the total stellar mass of a galaxy. Firstly, 
the rest-frame NIR is less sensitive to dust extinction than 
the UV and optical and therefore provides a more unbiased 
view of stars in galaxies. Secondly, the most massive evolved 
stars (i.e. those that persist on Gyr timescales) in galaxies 
are preferentially redder and emit strongly in NIR wave- 
lengths so the infrared light provides a better representa- 
tion of the high-mass end of the stellar mass function which 
will contribute significantly to the total stellar mass. The 
UV/optical wavebands are dominated by light from young 
luminous short-lived stars which contribute significantly less 
to the total stellar mass. Finally, and particularly relevant 
for our sample, in spectroscopic surveys covering a wide 
redshift range, the availability of additional photometry at 
longer wavelengths allows us to sample the same rest-frame 
portion of the galaxy SED at high redshifts as sampled by 
the UV and optical filters at low redshifts. For these reasons, 
rest-frame NIR mass-to-light ratios have o ften been used to 
derive robust stellar masses for galaxies (|Bell et all I2OO3I : 
iDrorv eFa l. 2004). Although it is true that stellar popula- 
tion synthesis (SPS) models suffer from larger uncertainties 
in modelling the rest-frame NIR SEDs of galaxies, it is im- 
portant to quantify the effect these uncertainties have on 
the stellar masses and how these trade off with the inclusion 
of additional data points in the fitting procedure. 

We use the FAST code to derive stellar masses 
for a) the 27,305 WiggleZ galaxies unmatched to wide- 
field IR surveys (FUV, NUV,ugriz photometry), b) the 
11,919 WiggleZ galaxies at 0.3<z<1.0 that are also de- 
tected at >5(7 in at least one of the UKIDSS bands 
{FUV,NUV,ugrizYJHK photometry) and c) the 6117 



galaxies in the same redshift range that are also detected 
in WISE {FUV,NUV,ugrizYJHK,Wl,W2 photometry). 
About a quarter of the 11,919 UKIDSS matched galaxies 
are detected in all four UKIDSS bands - Y, J, H, K. In some 
cases, data is missing from one or more of the UKIDSS bands 
due to the fact that there are regions of sky that haven't 
been imaged in all four bands. However, ~20% of the galax- 
ies are only detected in the Y-band on account of being very 
blue and fall below the S/N threshold used to construct the 
UKIDSS catalogues in the redder UKIDSS bands. For this 
subset, it is also interesting to assess the effect of using the 
low S/N fluxes in the red bands in the SED flts versus ig- 
noring the galaxy photometry in these bands. 

In all cases we assume BC03 SPS models, an exponen- 
tially decaying star formation history (SFH) with r ranging 
from 0.01 to 30Gyr - the latter representing constant star 
f ormation - a grid of metallicities with 0.001<Z<0.05 and 
a lChabried (|2003l ) IMF. Dust extinction is allowed to vary 
over 0<Av <2. We will test how some of these parameter 
choices can affect the stellar mass estimates, later in this 
section. 

4.1.1 Differences in IR undetected, UKIDSS detected and 
UKIDSS+WISE detected sub-samples 

In Figure [3l we plot the stellar mass distributions for the IR- 
missed, UKIDSS and UKIDSS-H WISE samples as a function 
of the FUV absolute magnitude. The FUV absolute magni- 
tude used throughout this paper is derived by k-correcting 
the observed NUV magnitude for each galaxy using a Ly- 
man Break Galaxy (LBG) template reddened by Av=0.14 
mags. This LBG template produces a very good match to 
the aver age observed (NU V — r) colours of the WiggleZ 
galaxies (|Blake et al.ll200^ ). Some of the WiggleZ galaxies, 
particularly those detected in the IR, are redder than this 
template, but our best-fit SEDs only begin to diverge from 
the LBG template at rest-frame wavelengths of A >2000A. 
The FUV luminosity is found to be relatively insensitive to 
the choice of template and we find that differences in the 
FUV luminosity computed using this LBG template ver- 
sus the individual best-fit SEDs, are at most ~10%. As we 
will be deriving stellar masses and SED fit parameters us- 
ing a range of different inputs, when considering how these 
fit parameters change as a function of rest-frame galax;y lu- 
minosity or colour, it is important to ensure that the rest- 
frame quantity used does not depend critically on the form 
of the best-fit SED. For this reason, the FUV luminosity 
is used throughout this paper. Also, for ease of computa- 
tion, the single LBG template is used to derive the FUV 
k-corrections and we have checked that using the individual 
SEDs instead does not affect any of our conclusions. 

We find that the IR detected WiggleZ galaxies have 
a fairly tight distribution in stellar mass while the less lu- 
minous galaxies that are unmatched to wide-field IR sur- 
veys, span a much larger range in stellar mass and have 
in general lower stellar masses and larger errors on the in- 
dividual stellar mass estimates. The median stellar masses 
are logio(M./M0)=9.6±O.7, 10.2±0.5 and 10.4±0.4 for the 
IR missed, UKIDSS matched and UKIDSS+ WISE matched 
sub-samples respectively. It is important to clarify the mean- 
ing of the various errors on the stellar masses quoted 
throughout this paper. The quoted errors on the median 
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masses, refer to the standard deviation for the entire sam- 
ple. The median la errors on the individual mass estimates, 
derived using Monte Carlo simulations, is typically slightly 
smaller - ~0.3-0.4 dex. The error on deriving this median on 
the other hand, assuming Gaussian statistics, is extremely 
small given our large sample size (<0.01 dex). From hereon, 
the errors on the median reflect the standard deviation of 
the sample, but we will also refer to the la errors computed 
from the Monte Carlo simulations as a means of assessing 
the accuracy of the stellar mass estimates. 

We also note that the IR undetected galaxies seem to 
have a bimodal distribution in the stellar mass - Mfuv plane 
with a cloud of galaxies with similar stellar masses and FUV 
luminosities to the IR-detected galaxies, and a separate 
cloud of fainter, lower stellar mass galaxies. There is a clear 
separation between these two populations at a stellar mass 
of logio(M,/M0)~ 9.7. The lower stellar mass galaxies are 
predominantly at lower redshifts whereas at higher redshifts, 
the WiggleZ colour cuts seem to select more massive galax- 
ies. We note however that the IR detected population, which 
is predominantly more massive than logio(M*/M0)~ 9.7, 
contains galaxies over the entire redshift range of the Wig- 
gleZ survey. Using the IR detected galaxies, we check that 
excluding the IR data from the SED fitting, does not repro- 
duce this observed bimodality. We find therefore that at the 
high redshift end of the WiggleZ sample, there is a signifi- 
cant population of massive galaxies that are undetected in 
wide-field NIR surveys, whereas at the low redshift end of 
the sample, those galaxies that are undetected in the NIR, 
are predominantly fainter and have lower stellar masses. We 
note further that this bimodality is also somewhat mirrored 
in the observed colour-magnitude diagram in Figure[2]where 
we see that some of the WiggleZ galaxies that are undetected 
in the NIR surveys, are just as red in terms of their optical 
colours as the NIR detected galaxies. 

Reliable age constraints are difficult to obtain from 
these SED fits. For single component star formation his- 
tories used in the FAST code, the age refers to the age of 
the galaxy since the onset of star formation, and ranges from 
~250 Myr for the IR-undetected population, with typical la 
errors of 2-3 Gyr on these age estimates, versus ~400-1000 
Myr ±1.5 Gyr for the IR detected galaxies. The age is very 
degenerate with the dust extinction with both parameters 
having a similar effect on the observed colours of galaxies. 
The UKIDSS and WISE detected galaxies are found to have 
higher Ay by ~ 0.4 magnitudes, compared to the galaxies 
not detected in the infrared. However, the typical la errors 
on these Ay estimates are ~0.6 magnitudes so the age and 
dust extinctions are not particularly well constrained. 
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Figure 3. Stellar mass as a function of FUV absolute magni- 
tude for the WiggleZ galaxies not detected in wide-field IR sur- 
veys (top), those detected in UKIDSS out to 2.2/im (middle), and 
those detected in both UKIDSS and WISE out to 4.6Atm (bot- 
tom). The grcyscalc represents the density of points in each panel. 
The galaxies detected in the IR have a much tighter stellar mass 
distribution compared to the bluer galaxies missed in wide-field 
IR surveys. They are also clearly more massive. In each panel, we 
show the best-fit straight line through these points as the solid 
line to guide the eye and the best-fit lines for the other panels 
are shown as the dashed lines for comparison. For the top panel 
of galaxies undetected in the IR, we have fit two separate best-fit 
lines to the two clouds of galaxies at logio(M«/MQ)< 9.7 (blue 
line) and logio(M,/M0)> 9.7 (cyan line). 



4- 1-2 How does the IR data change the stellar masses? 

Are these discrepancies between the stellar masses of the IR 
undetected and detected population, due to differences in 
the physical properties of these galaxies? Or do they simply 
result from including different photometric bands in the fits 
for the two samples? In order to assess this, we compute 
masses for the 11,919 galaxies matched to UKIDSS, but ex- 
cluding the NIR data from the SED fitting - i.e. we want to 
look at the effect of using only the UV/optical photometry 
in the SED fits for the exact same galaxies. The la errors on 
the stellar mass estimates, derived using FAST, are plotted 



as a function of the stellar mass in the bottom two panels of 
Figure |4] both when including and excluding the NIR data 
from the SED fits for the subset of 11,919 galaxies. 

Figure|4]clearly demonstrates that the stellar masses are 
better constrained with considerably smaller la errors when 
we include the NIR data in the SED fits in the FAST code. 
For the 11,919 UKIDSS detected galaxies, the median stellar 
mass is logio(M,/M0) = lO.2±O.5 when the UKIDSS data is 
included in the fitting versus logio(M»/M0) = lO.4±O.6 when 
the UKIDSS data is excluded. Once again we note that these 
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errors refer to the total standard deviation for the sample 
and that the error on the medians are considerably smaller, 
making these differences in the two mass estimates highly 
statistically significant. As the removal of NIR photometry 
from the SED fits only serves to make the median mass 
of the NIR detected population even larger, we conclude 
that the difference in mass seen between the NIR detected 
and undetected galaxies in Section 14.1. II is real and not 
an artefact of including different bands in the SED fitting. 
The NIR detected WiggleZ galaxies are therefore on average 
more massive than those not detected in the infrared. 

We now turn our attention to considering whether the 
addition of the NIR has indeed led to better constraints 
on our stellar mas ses as suggested by Figure [l] Recently, 
iTavlor et al.l (|201ll ) have argued that including near infrared 
data from the UKIDSS-LAS when SED fitting to galaxies in 
the GAMA survey, results in stellar mass estimates that are 
highly discrepant between the optical and optical+NIR fit 
cases. They find that the optical+NIR derived masses are 
inconsistent with the optical only masses at the >3cr level for 
~25% of the GAMA galaxies. They conclude that the NIR 
can therefore not be used to provide reliable stellar mass 
constraints for their dataset. We conduct a similar test in 
the case of our WiggleZ galaxies using the FAST outputs for 
our stellar masses. Taking into account the formal Icr errors 
on the individual mass estimates derived using the Monte 
Carlo simulations, we look at the fraction of galaxies where 
the masses from the optical and optical+NIR runs, are con- 
sistent within these Icr errors. We find that this is true in 
>70% of the galaxies indicating that although there is a 
systematic shift in the median mass towards lower masses 
when including the NIR data, the masses thus obtained are 
still consistent with the larger errorbars derived using the 
UV and optical data only. The fraction of galaxies with 
>3(T discrepant mass estimates between the optical and op- 
tical+NIR fits, is only ~4% in the c ase of the FAST out puts 
compared to the ~25% reported in lTavlo~ et all (|201ll ) for 
the GAMA galaxies. We therefore find that the NIR data 
does help better constrain the SED models for the majority 
of the WiggleZ galaxies. 

Can the differen ces between our conclusions and those 
of lTavlor et~all (|201ll ') , be explained by FAST's use of a tem- 
plate error function to model uncertainties in stellar popula- 
tion synthesis models over certain wavelength ranges? This 
template error function is lowest in the rest-frame optical 
and increases both at rest-UV wavelengths due to the effects 
of varying dust extinction, and in the NIR due to uncertain- 
ties in the stellar isochrones (|Marastonl |2005| ). We recom- 
pute stellar masses for the UKIDSS detected galaxies with 
and without the near infrared data, but this time turning 
off the template error function within FAST. The template 
error functio n increases the phot ometric uncertainties at 
A > 10000 A (|Brammer et al.ll2008l ) and therefore the resuff- 
ing errors on the masses. Once again, we find that the stellar 
masses are more tightly constrained with the inclusion of 
the NIR data. The fraction of galaxies with >3o- discrepant 
mass estimates between the optical and optical+NIR fits 
now goes up to ~6%, still consider ably smaller than the 
25% reported bv lTavlor et~al] (|201ll ). Our improved stellar 
mass constraints with the NIR data, are primarily driven 
by the fact that the Icr upper bound on the stellar mass 
is brought down by 0.4-0.5 dex with the inclusion of the 



NIR photometry. The lower bound on the stellar mass com- 
puted using the Monte Carlo simulations, however remains 
unchanged. The NIR photometry therefore helps to rule out 
very large stellar masses for the majority of WiggleZ galax- 
ies and thereby tightens the stellar mass estimates. We also 
check that using fewer Monte Carlo simulations to calibrate 
the error estimates, leads to larger discrepancies (~6% of 
galaxies with inconsistent masses at the >3cr level) between 
the optical and optical+NIR mass estimates. This is to be 
expected given that as we decrease the number of Monte 
Carlo simulations, the errors also become less reliable. 

We check our results using the KG04 code used in con- 
junction with the PEGASE.2 SPS models. The KG04 code 
makes different assumptions regarding the priors on the in- 
put parameters for the SED fits and a key motivation of the 
current work is to test the robustness of stellar mass esti- 
mates to these differences in SED fitting codes. In contrast 
to the results obtained using FAST, we find little difference 
in both the median stellar mass estimate and the Icr error 
on it with and without the addition of the NIR photome- 
try. We use the KG04 code with both single component star 
formation histories and the addition of secondary bursts on 
top of the smooth underlying SFH. In both cases, the me- 
dian stellar mass remains unchanged on addition of the NIR 
data, as does the typical Icr error. More than 80% of galax- 
ies now have stellar masses that are consistent between the 
optical and optical+NIR SED fits within the Icr errors. Less 
than 2% have stellar masses that are discrepant at the >3cr 
level. 

Our study has therefore illustrated two key points. The 
effect that the NIR data has in constraining the stellar 
masses of galaxies, depends on the assumptions made in the 
model being fit. Hence, we get slightly different results with 
the FAST and the KG04 code although it is important to 
stress that with the NIR photometry included, the median 
masses produced by the two codes are very similar. Secondly, 
it is important to have well calibrated errors on the stellar 
masses that take into account both the error on the photom- 
etry as well as the error in the templates in the NIR. If we 
do not include the template error function within FAST the 
discrepancy between the optical and optical+NIR computed 
masses increases. Similarly, we found that if we use a smaller 
number of Monte Carlo simulations to calibrate the errors, 
the discrepancy once again increases. We conclude that the 
NIR data can for some models lead to tighter constraints 
on the stellar masses, and find no evidence for widely dis- 
crepant stellar masses computed using the optical only and 
optical+NIR data with either the FAST or KG04 code. 

4-1.3 Effect of adding in low signal-to-noise fl,uxes 

There are a significant number of blue galaxies in WiggleZ 
that are detected in the UKIDSS y-band but are faint and 
at low S/N in the redder UKIDSS bands. In these cases, 
we have so far ignored the JHK photometry from UKIDSS 
in the SED fitting. We now assess whether including these 
low S/N fiux measurements results in any improvement in 
the SED fits and any difference to the inferred fit parame- 
ters. In order to do so, we select a subset of ---^1000 galaxies 
that are detected in the K-band but undetected in H and 
K. We perform forced photometry on these galaxies in the 
H and /('-bands using fixed apertures centred on the Wig- 
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Figure 4. The distribution of la error on the stellar mass esti- 
mates as a function of stellar mass, computed using FAST, for the 
sample of 11,919 WiggleZ galaxies matched to UKIDSS when us- 
ing different photometric bands in the SED fitting. The greyscale 
represents the density of points in each panel while the points 
with error bars represent the median and standard deviations of 
the samples. The Icr errors are systematically larger when the NIR 
data is excluded from the fits (bottom panel) confirming that the 
inclusion of this data where available, provides better constraints 
on the stellar masses (middle panel). The stellar masses too are 
over-estimated when the NIR data is excluded from the fits. The 
GALEX UV-bands however make little difference to the inferred 
median stellar mass and errors (top and middle panels) although 
there is evidence for a cloud of galaxies where the stellar masses 
are more poorly constrained without the UV data. This cloud of 
galaxies is seen to migrate from the top left of the top panel to 
the bottom right of the middle panel in the figure. 



gleZ positions. The fluxes and errors are measured in several 
difTerent apertures and using the high S/N detections, we 
find that a 4" aperture provides a reasonable match to the 
UKIDSS catalogue Petrosian magnitudes of the galaxies. At 
low signal-to-noise however, a larger aperture often results 
in more noisy photometry so we choose the aperture flux 
estimate with the highest signal-to-noise ratio to represent 
the total flux from the galaxy. These aperture sizes range 
from 2 to 5". The median S/N is ~3 in both the H and K- 
bands. Stellar masses are then calculated using FAST but 
now including the results from the forced photometry in the 
H and K bands. In ~60% of the galaxies, we find an im- 
provement in the Icr stellar mass errors by a factor of ~2, 
when including the forced photometry fiuxes. However, the 
inferred stellar masses themselves change only by ~0.01 dex. 
In ~5% of the galaxies, the inclusion of the forced photome- 
try fluxes considerably worsens the SED fits. These galaxies 
are, as expected, typically < Icr detections in the H and 
if-band images, and the inclusion of these very noisy fluxes 
worsens the quality of the fits. We conclude that using low 
signal-to-noise aperture fluxes versus ignoring the photom- 
etry in these bands, does not materially change the stellar 
mEiss estimates for the majority of the galaxies although it 
can help to reduce the errors on the individual stellar mass 
estimates. 

4-1-4 How does the UV data change the stellar masses? 

We also consider the effect of removing the GALEX UV pho- 
tometry from the SED fits for the 11,919 galaxies matched 
to UKIDSS. The corresponding stellar masses and errors 
computed using FAST and the BC03 models, are plotted in 
the top panel of Figure [3] We find that the median stellar 
mass goes down by only ~0.02 dex versus when the GALEX 
photometry is included (middle panel of Figure [Jjl , and the 
median error on this stellar mass remains unchanged. More 
than 80% of the galaxies have stellar masses that are consis- 
tent within the Icr errors when considering the outputs with 
and without the UV photometry. 

However, we find in Figure |4] that there is a cloud of 
galaxies where the removal of the GALEX photometry ac- 
tually leads to larger errors on the stellar mass estimates. 
We investigate this cloud of galaxies in more detail. Without 
the UV photometry, the best-fit SEDs suggest that these are 
young, highly star-forming galaxies with significant amounts 
of dust extinction. However, when the UV photometry is in- 
cluded, the best-fit SEDs suggest that these galaxies are on 
average more massive, significantly older and less dusty and 
with significantly lower star formation rates than inferred 
from SED fits without the UV data. In Figure|4l these galax- 
ies migrate from the top left portion of the top panel, to- 
wards the bottom right portion of the middle panel. The UV 
photometry is dominated by light from instantaneous star 
formation rather than reflecting the total underlying stellar 
mass. In most galaxies, it therefore provides little constraint 
on the total stellar mass. However, we see that for some of 
the WiggleZ galaxies, the lack of UV photometry can lead 
to an over-estimate of the dust extinction and star forma- 
tion rate and an under-estimate of the age of the galaxy. 
As a consequence of the SEDs being flt by younger stellar 
populations, the total stellar mass is also under-estimated 
without the UV data, in these galaxies. 
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We also consider the efTect of the UV photometry on the 
27,305 fainter, bluer galaxies that are undetected in the IR 
and therefore where NIR photometry is no longer included 
in the SED fitting. Once again, there is little effect on the 
median stellar mass, which goes down by O.Of dex on re- 
moving the GALEX bands from the fits, and the median la 
errors remain unchanged. Almost 90% of galaxies now have 
stellar masses that are consistent within the Icr errors when 
considering the outputs with and without the UV photome- 
try. These 27,305 galaxies have no NIR data and in the case 
of these galaxies, we find no evidence for a separate cloud of 
galaxies where the lack of UV photometry leads to poorer 
constraints. 

4.2 Effect of Stellar Population Synthesis Models 

We now consider the effect of different SPS models on the 
stellar mass estimates with all other parameters in the SED 
fitting held constant. We use a Salpeter IMF and a set of sim- 
ple single component exponentially decaying SFHs. FAST 
allows fitting of both the BC03 and iMarastonI (|2005l ') mod- 
els to the multi-wavelength data. While the latter includes 
a prescription for TP-AGB stars, the BC03 models do not. 
The contribution of TP-AGB stars to the NIR colours of 
galaxies is still widely debated with local galaxies showing 
significant dependen ces of the colours on this component 
minian et al.ll200^ ) while the contribution in higher red- 
shift galaxies is deemed to be less significant (|Kriek et al.l 
I2OIOI '). TP-AGB stars are only expected to affect the NIR 
colours of intermediate age (~lGyr) stellar populations. We 
also derive stellar masses using the KG04 code and a set of 
PEGASE.2 SEDs, in which nebular emission lines can be 
turned on or off in order to consider their effect on the stel- 
lar masses. Note that although the KG04 code incorporates 
the facility to include secondary bursts of star formation, 
for this particular test, only single component SFHs are al- 
lowed so as to allow direct comparison with the results from 
FAST where secondary bursts cannot be included. Having 
already demonstrated that the NIR data can help in better 
constraining the stellar masses for FAST, and that it does 
not materially change the stellar masses for the KG04 code, 
we use the UV, optical and NIR photometry in the SED fits 
in all cases where the NIR is available. We keep the SFH 
and IMF fixed so the effect of the different models on the 
stellar mass estimates can be isolated. 

For the NIR detected population, we find the reduced 
values and Icr errors on the stellar masses to be very sim- 
ilar for the BC03 and Maraston models. However, the me- 
dian stellar mass is ~0.2 dex lower when using the Maraston 
models. For the bluer subset of galaxies undetected in the 
NIR, the average difference between the Maraston and BC03 
derived stellar masses now falls to ~0.1 dex. Figure [5] shows 
the distribution of stellar mass as a function of }Afuv for 
both WiggleZ galaxies detected and undetected in the NIR 
for each of the different SPS models studied. Changing from 
the BC03 models to the Maraston models primarily affects 
the red cloud of galaxies detected in the IR. This is expected 
given that we infer older ages for these IR detected galaxies 
suggesting that they have slightly more evolved stellar pop- 
ulations where TP-AGB stars may become more important. 

Before considering the differences between the 
BC03/Maraston and PEGASE.2 models, we first assess the 



effect that nebular emission lines have on the stellar masses, 
by computing these masses with the nebular emission 
switched on and off in the PEGASE models. We find that 
turning off the nebular emission in the models increases 
the stellar masses as expected by raising the best-fit stellar 
continuum. However, the median increase in the stellar 
masses due to lack of nebular emission in the models, is only 
~0. 01/0. 02 dex for those galaxies detected in the IR and 
those that are undetected respectively. The WiggleZ galax- 
ies are selected to have strong emission li nes in order to aid 
the a cquisition of redshifts for cosmology jDrinkwater et all 
|2010D . However, even in these strong emitters, the impact 
of including nebular emission in the models on the stellar 
masses, appears to be negligible for the majority of the 
galaxies and smaller t han the offsets typically seen in 
high redshift LBGs (e.g. Ide Barros et al.l [20121 : IStark et all 
l2012h . We note however, that the prescription for nebular 
emission in the PEGASE.2 SEDs, is relatively simple and 
does not for example allow for variations in the line ratios. 
There are of course WiggleZ galaxies where the emission 
lines make a significant contribution to the broadband flux, 
and where the inclusion of nebular emission in the models, 
results in a more significant offset to the stellar mass. Some 
example SED fits of galaxies with strong nebular emission 
can be found in Appendix A. However, on average, the 
effect of including emission lines in the models, on the 
stellar masses, is found to be very small. 

Having demonstrated that the inclusion of nebular 
emission lines in the SPS models, appears to have little im- 
pact on the median stellar masses of the WiggleZ galax- 
ies, we now look at the differences in stellar mass between 
the PEGASE.2 models and the Maraston and BC03 mod- 
els. The results can be seen in Figure [5] We find that the 
PEGASE.2 models always produce stellar masses that are 
higher than those derived by the BC03 and Maraston mod- 
els. These masses are higher by ~0.1 (0.2) dex with respect 
to BC03 (Maraston) for the galaxies detected in the IR, 
and ~0.2 (0.3) dex with respect to BC03 (Maraston) for the 
bluer galaxies undetected in the IR. We have ruled out neb- 
ular emission as the cause for these mass differences, and 
shown that the presence of emission lines would decrease 
rather than increase the stellar masses computed with the 
PEGASE models. These differences, which seem to be ac- 
centuated at faint FUV luminosities, must therefore result 
from differences in the input stellar libraries in each of the 
models and the fact that these stellar libraries traverse dif- 
ferent regions in physical parameter space. We note that 
restricting the PEGASE models to solar metallicity only, 
produces slightly more consistent stellar masses between PE- 
GASE and BC03 by reducing the PEGASE produced steUar 
masses by ~0.1 dex. 

Unlike the stellar masses which are reasonably well con- 
strained, the best-fit ages are once again more poorly con- 
strained and show significant deviations between the models. 
Nebular emission is found to produce younger ages and the 
PEGASE.2 models generally produce higher median ages 
than BC03 and Maraston. However, we emphasise that the 
uncertainties on these age estimates are considerable. 



© 0000 RAS, MNRAS 000, 000-000 



12 




12p 
11.5 



NIR Detected 
NIR Undetected i 





Figure 5. Stellar mass versus Wlpuy for the WiggleZ galaxies undetected in the IR (blue dashed) and those detected in the IR (red 
solid) for three different choices of SPS models. Contours represent the density of points computed using a kernel density estimator, 
over a grid traversing the parameter space shown. The stellar masses are lowest for the Maraston models which include TP-AGB stars, 
followed by BC03 and then PEGASE, which produces the largest masses. The discrepancy between PEGASE and BC03/Maraston, is 
largest at faint Wlpuy and particularly noticeable for the blue cloud of IR undetected galaxies. The discrepancy between BC03 and 
Maraston on the other hand affects more the IR detected, bright yipjjv galaxies that populate the top left corner of the panels. 



4.3 Effect of Different Initial Mass Functions 

The form of the stellar initial mass function assumed dur- 
ing SED fitting, can have a significant effect on the inferred 
stellar masses of galaxies. In this section, we systematically 
quantify the difference in stellar masses using different IMFs 
but assuming the same SPS models, photometric filters and 
star formation histories. We use the subset of 11,919 Wig- 
gleZ galaxies matched to UKIDSS, for the comparison. We 
change the IMF choice in both FAST and KG04 while keep- 
ing the rest of the input parameters fixed, in order to quan- 
tify the resulting change in the median stellar masses of the 
WiggleZ galaxies. 

Firstly, using FAST with the BC03 models and expo- 
nentially decaying SFHs, we find that the median difference 
between the Salpeter and Chabrier IMFs is ~0.24 dex. The 
Chabrier IMF has the same power-law slope as Salpeter at 
the high-mass end but turns over at M<1Mq. 

FAST is also used with the Maraston models and expo- 
nentially decaying SFHs t o look at differe nces in stellar mass 
between the Salpeter and lKroupal (|200lh IMF. The Kroupa 
IMF predicts stellar masses that are on average ~0.20 dex 
lower than Salpeter. It too has the same power-law slope as 
Salpeter at the high-mass end but turns over at M< O.8M0. 
These differences in stellar mass due to the IMF, show little 
dependence on the FUV luminosity or galaxy colour. 

Finally, we use the KG04 code with PEGASE mod- 
els and multi-component SFHs, to look at the difference 
between the Salpeter, Kroupa, BG03 and Chabrier IMFs 
and to check if these differences are enhanced or reduced 
when additional bursts of star formation are allowed. The 
BG03 IMF has a shallower slope than Salpeter at the high- 
mass end and turns over at M<O.5M0. It, results in stellar 
masses that are on ave rage 0.24 dex lower than Salpeter 
l|Glazebrook et al.|[20o3 ). The stellar masses obtained with 
the BG03 and Chabrier IMFs are remarkably similar with a 
median difference of <0.001 dex between them. The Kroupa 
IMF results in stellar masses that are ~0.1 dex higher than 
those produced by BG03 and Chabrier. 

The simple power-law form of the Salpeter IMF is 
known to overestimate the number of l ow-mass stars in 
galaxies and therefore t he stellar mass l|Cappellari et all 
l2006l : lFerreras" et al.|[2"008h and more recent analytical forms 



of the IMF such as those considered in this work, have, by 
design a more realistic break in the IMF at low masses and 
therefore predict lower stellar masses on average. We find 
that the Kroupa IMF results in stellar masses that are ~0.1- 
0.2 dex lower than Salpeter while the BG03 and Chabrier 
IMFs produce very similar stellar masses that are ~0.24 dex 
lower than Salpeter and up to 0.1 dex lower than Kroupa. 



4.4 Effect of Different Star Formation Histories 

4-4- 1 Single Component Star Formation Histories 

In order to study the dependence of the stellar mass esti- 
mates on the assumed SFH, we begin by using FAST to fit 
BC03 models with three different star formation histories to 
the photometric data. These correspond to a simple expo- 
nentially decaying SFH, a truncated SFH which corresponds 
to constant star formation between tonset and tonsct+i", and a 
delayed exponentially decaying star formation history with 
SFR~ t exp(-t/r). We note that the FAST code used for 
the stellar mass estimates currently does not incorporate 
the ability to include exponentially increasing or inverted-r 
models which have recently been argued to provide a bet- 
ter description of the extinction and star formation rates 
inferred from rest-fram e UV data alone for galaxies at z>l 
(|Maraston et al.ll201ol ). 

In Figure[n]we show the distribution of stellar masses for 
the 11,919 UKIDSS detected WiggleZ galaxies as a function 
of Mfuv for each of the three SFHs. We find that with a 
large sample of spectroscopically confirmed galaxies such as 
ours, the stellar mass is insensitive to the choice of star for- 
mation history for simple single component SFHs. However, 
as might be expected for these very blue emission-line galax- 
ies, the truncated SFH produces relatively poor reduced-x^ 
values as it is typically invoked to describe the SFHs of more 
pas sive systems such as quiescent and p ost-starburst galax- 
ies (jPaddi et al.ll2004l : iKriek et al.ll2010l ). 

Although the stellar masses are insensitive to the choice 
of SFH, the ages and dust extinctions, which are in general 
more poorly constrained, show significant variations as the 
SFH is changed and are not particularly well constrained. 
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Figure 6. Stellar mass distribution as a function of yipuv ^ for 
the IR detected WiggleZ galaxies, computed with various different 
prescriptions for single component SFHs. The contours represent 
the density of points estimated using a kernel density estimator 
over a grid covering the parameter space shown. The stellar mass 
is seen to be insensitive to the choice of star formation history for 
single component SFHs. 



Figure 8. Difference in stellar mass obtained using the KG04 
code and PEGASE models on inclusion of bursts, versus the burst 
mass fraction for all 39,701 WiggleZ galaxies. The greyscale de- 
notes the density of points. The median burst fraction is ~1 per 
cent and the inclusion of the bursts in the SED models, adds 
~0. 1-0.3 dex of stellar mass on average to the galaxies. 



4. 4-2 Addition of Bursts 

Although the stellar mass is found to be insensitive to the 
assumed star formation history for a single component SFH, 
it is interesting to quantify the effect of more complex SFHs 
on the mass estimates. FAST does not allow for addition of 
random bursts of star formation on top of the single com- 
ponent star formation history. In order to assess the effect 
of these bursts on the stellar masses, we use the KG04 code 
with the PEGASE. 2 SEDs including nebular emission and 
assuming a BG03 IMF. As demonstrated in Section [4.31 the 
choice of IMF has little effect on the stellar masses for IMFs 
with realistic breaks at low masses. We examine the impact 
of starbursts on the stellar masses, both for the IR detected 
subset of 11,919 galaxies as well as the bluer 27,305 galaxies 
that are undetected in wide-field IR surveys. 

The results are shown in FigureQwhere we plot the stel- 
lar masses as a function of yipuv for both subsets of galax- 
ies with and without additional bursts of star formation. We 
find that the addition of bursts of star formation leads to 
an increase in the stellar masses in both cases. This increase 
is more pronounced for those galaxies not detected in the 
IR where the addition of bursts adds ~0.3 dex to the stellar 
masses, versus 0.1 dex for galaxies detected in the IR. For 
both subsamples, the increase in stellar mass on inclusion of 
bursts is also more pronounced at higher FUV luminosities. 
This is because the secondary bursts essentially "hide" the 
more massive evolved stellar populations by dominating the 
UV-light. 

We check these results using the MAGPHYS code, 
which makes different assumptions in the implementation 
of the secondary burst as discussed in detail in Section O 
As previously stated, MAGPHYS is not readily apphca- 
ble to the large sample sizes assembled here, but, unlike 
FAST, it does include additional bursts of star formation. 
We therefore compare the FAST and MAGPHYS derived 
stellar masses for a randomly selected sub-sample of 800 



WiggleZ galaxies that are also detected in UKIDSS and 
a similar sub-sample of WiggleZ galaxies not detected in 
UKIDSS. We find that the mean difference in stellar mass 
between FAST and MAGPHYS is 0.1 dex for the UKIDSS 
detected galaxies and 0.25 dex for the UKIDSS undetected 
galaxies, consistent with the results derived using the KG04 
code due to the effect of bursts. 

In Figure |8] we plot the difference in stellar mass on in- 
clusion of bursts versus the burst mass fraction derived using 
KG04 for all ~40,000 WiggleZ galaxies. Most of the galaxies 
have typical burst fractions of ~1% - i.e. the mass produced 
in the secondary burst is just 1% of that produced through 
continuous star formation. However Figure |8] clearly shows 
that the difference in stellar mass between the burst and no 
burst SED fits, is not particularly well correlated with the 
burst fraction. This demonstrates that allowing secondary 
bursts in the SED fits, does not just change the total stellar 
mass of the galaxy but also results in best-fit SEDs with 
different properties - age, Ay, t - compared to when bursts 
are not included in the SED fits. 

For multi-component SFHs, the ages output are the 
mass-weighted ages of the different stellar populations so a 
direct comparison to the ages produced by fitting single com- 
ponent SFHs, cannot be made. The median mass weighted 
ages are ~2-3 Gyr. The addition of secondary bursts there- 
fore produces older galaxies than when these bursts are not 
included. This is because the underlying primary stellar pop- 
ulation being hidden by the current burst is also older than 
when the bursts are not included. 

We also consider the variation of the best-fit dust ex- 
tinction parameter Av with FUV luminosity and stellar 
mass. We consider the outputs from the KG04 code includ- 
ing both nebular emission as well as secondary bursts of star 
formation. Although these dust extinctions are typically not 
very well constrained with median la errors of ~0.2 mags, 
we can nevertheless assess whether there is any evidence for 
a trend in best fit Av with the FUV luminosity or stel- 
lar mass. Figure |9] shows the median best-fit Ay in bins 
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Figure 7. Stellar mass distribution as a function of FUV absolute magnitude for WiggleZ galaxies detected and undetected in infrared 
surveys. The greyscale denotes the density of points while the individual points with error bars represent the median and standard 
deviations of the samples. Stellar masses are calculated using the KG04 code, PEGASE SEDs, assuming a BG03 IMF and with additional 
bursts of star formation turned on or off. Additional bursts of star formation add ~0.1 dex in stellar mass to galaxies detected in the 
NIR and ~0.3 dex of stellar mass in galaxies undetected in the NIR. The effect is more pronounced at brighter FUV luminosities. 



of Lfuv and stellar mass for all ~40,000 WiggleZ galaxies 
with the shaded regions representing the la errors from the 
SED fits. Despite the large errors, we find a strong trend 
with stellar mass with the most massive galaxies also be- 
ing fit to have larger dust extinction values. The trend with 
FUV luminosity is less strong but there is some evidence 
for WiggleZ galaxies with lower FUV luminosities also hav- 
ing higher best-fit Av- These trends are confirmed from the 
outputs derived using the FAST code which includes only 
single component SFHs and no nebular em ission lines in the 
SPS models. Recently, iBuat et al.1 (|2012h observed a simi- 
lar trend in dust attenuation using Herschel data to analyse 
a sa mple of intermediat e redshift galaxies in the GOODS 
field. iHeinis et al.l (|2012l ) too find a slight increase in dust 
attenuation at low FUV luminosities although the trend is 
essentially fiat at Lfuv > 1O^''L0, once again consistent 
with our data in Figure |9l However, given the biased colour 
selection of WiggleZ targets and the considerable uncertain- 
ties and degeneracies involved in accurately constraining Av 
from SED fits to the available broadband photometry for 
this sample, we caution against over interpretation of these 
observed trends. 



4.5 Summary & if-Band Mass-to-Light Ratios 

We have quantified the sensitivity of stellar mass estimates 
for the very blue population of WiggleZ galaxies at 0.3 < 
2 < 1.0 to changes in various input parameters and assump- 
tions made during the SED fitting. In particular, we have 
looked at how the median stellar masses of this large sample 
are affected by changes to the input photometry, choice of 
SPS model, choice of IMF and choice of star formation his- 
tory. The SPS models and the addition of secondary bursts 
result in the biggest changes to the median stellar masses - 
of the order of 0.3 dex for these very blue galaxies. Table [T] 
summarises the results of the various SED fits carried out 
in this Section and provides median stellar masses for dif- 
ferent subsamples of the WiggleZ population output using 
different codes and different input parameters to the SED 
fitting. A key result of this paper is that the stellar masses 
seem to be relatively robust to changes in the various input 
parameters to the SED fitting as well as the choice of SED 
fitting code. 
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Figure 9. Median best-fit dust extinction parameter Ay derived using the KG04 code with the PEGASE.2 SPS models, as a function of 
FUV luminosity and stellar mass for all ~40,000 WiggleZ galaxies. The points represent the median values in bins of hpuy and stellar 
mass with the shaded regions denoting the la errors on from the SED fitting. Despite the large errors, we note a trend of increasing 
dust extinction with increasing stellar mass and decreasing ^fuv ™ this sample. 
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Table 1. Summary of median stellar masses of WiggleZ galaxies computed using different SED fitting codes and when changing different input parameters in the SED fitting. 



Code 


Ngal 


Filters 


SPS Model 


Nebular Lines 


SFH 


IMF 


<logio(M./Mo)> 


FAST 


27,305 


GALEX,ugriz 


BC03 


No 


Exp decaying 


Chabrier 


9.6±0.7 


FAST 


27,305 


ugriz 


BC03 


No 


Exp decaying 


Chabrier 


9.5±0.6 


FAST 


11,919 


GALEX,ugrizYJHK 


BC03 


No 


Exp decaying 


Chabrier 


10.2±0.5 


FAST 


11,919 


ugrizY J H K 


BC03 


No 


Exp decaying 


Chabrier 


10.2±0.4 


FAST 


11,919 


GALEX,ugriz 


BC03 


No 


Exp decaying 


Chabrier 


10.4±0.6 


FAST 


6117 


GALEX,ugrizYJHK,3A,4.6fj.m 


BC03 


No 


Exp decaying 


Chabrier 


10.4±0.5 


FAST 


11,919 


GALEX,ugrizYJHK 


BC03 


No 


Exp decaying 


Salpeter 


10.4±0.6 


FAST 


11,919 


GALEX,ugrizYJHK 


Maraston 


No 


Exp decaying 


Salpeter 


10.3±0.6 


KG04 


11,919 


GALEX,ugrizYJHK 


PEGASE 


Yes 


Exp decaying 


Salpeter 


10.5±0.4 


KG04 


11,919 


GALEX,ugriz 


PEGASE 


Yes 


Exp decaying 


Salpeter 


10.5±0.4 


KG04 


11,919 


GALEX,ugrizYJHK 


PEGASE 


No 


Exp decaying 


Salpeter 


10.5±0.4 


KG04 


27,305 


GALEX,ugriz 


PEGASE 


Yes 


Exp decaying 


Salpeter 


lO.OitO.6 


KG04 


27,305 


GALEX,ugriz 


PEGASE 


No 


Exp decaying 


Salpeter 


lO.OitO.6 


FAST 


27,305 


GALEX,ugriz 


BC03 


No 


Exp decaying 


Salpeter 


9.8±0.7 


FAST 


27,305 


GALEX,ugriz 


Maraston 


No 


Exp decaying 


Salpeter 


9.7±0.7 


KG04 


27,305 


GALEX,ugriz 


PEGASE 


No 


Exp decaying 


Salpeter 


lO.OitO.6 


FAST 


11,919 


GALEX,ugrizYJHK 


BC03 


No 


Delayed Exp 


Chabrier 


10.2±0.5 


FAST 


11,919 


GALEX,ugrizYJHK 


BC03 


No 


Truncated 


Chabrier 


10.2±0.5 


KG04 


11,919 


GALEX,ugrizYJHK 


PEGASE 


Yes 


Exp decaying + Burst 


BG03 


10.4±0.4 


KG04 


27,305 


GALEX,ugriz 


PEGASE 


Yes 


Exp decaying + Burst 


BG03 


10.0±0.7 


KG04 


11,919 


GALEX,ugrizYJHK 


PEGASE 


Yes 


Exp decaying + Burst 


Kroupa 


10.5±0.4 


KG04 


11,919 


GALEX,ugrizYJHK 


PEGASE 


Yes 


Exp decaying + Burst 


Chabrier 


10.4±0.4 


KG04 


11,919 


GALEX,ugrizYJHK 


PEGASE 


Yes 


Exp decaying + Burst 


Salpeter 


10.7±0.4 


FAST 


11,919 


GALEX,ugrizYJHK 


Maraston 


No 


Exp decaying 


Kroupa 


10.0±0.5 


MAGPHYS 


24 


GALEX,ugrizY J HK, 3 A,4:.6,12,22fj.m 


CB07 


No 


Exp decaying + Burst 


Chabrier 


10.8±0.3 


MAGPHYS 


800 


GALEX,ugrizYJHK 


CB07 


No 


Exp decaying + Burst 


Chabrier 


10.4±0.5 


MAGPHYS 


800 


GALEX,ugriz 


CB07 


No 


Exp decaying + Burst 


Chabrier 


9.9±0.7 
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We have also demonstrated that the near infra-red data 
can, for some codes and choice of priors, significantly reduce 
the 1(7 errors on the stellar masses by tightening the upper 
bound on the stellar mass estimate. We find no evidence 
that the near infra-red data actually worsens the quality 
of the SED fits for the WiggleZ galaxies with any of the 
codes. As previously mentioned, rest-frame near infra-red 
mass-to-light ratios have been used for some time to esti- 
mate the total stellar masses of galaxies. These mass-to-light 
ratios are significantly less dependent on galaxy colour than 
similar ratios derived at shorter wavelengths. Empirical re- 
lations of the dependence of stellar mass-to-light rati os on 
galaxy colour, have been derived by iBell et all (|2003h and 
are widely used in the literature. We therefore now com- 
pare this empirical relation in the near infrared K-hanA to 
the equivalent mass-to-light ratios derived from our SED 
fits for our sample of almost 12,000 WiggleZ galaxies that 
are detected in UKIDSS. Note that at the median redshift 
of our sample, the rest-frame isT-band is not sampled by 
the UKIDSS data. However, the availability of the UKIDSS 
photometry allows us to sample more of the galaxy SED 
than is sampled by the optical, making the extrapolation to 
the rest-frame K-hand more secure. We consider the out- 
puts from the KG04 code including secondary bursts of star 
formation and nebular emission lines in the SPS models. 



Age (Wlyr) 




g-r 



Figure 10. Stellar mass-to-light ratio in the rest-frame K-haad 
as a function of the observed (g — r) colour, derived from SED 
fits using the KG04 code and including both secondary bursts of 
star formation and nebular omission in the models. The points 
arc colour-coded according to the best-fit age of the galaxy. The 
dashed line shows the empir ical colour-based estimator of the 
M/Lx from lBell et al.l | |2003|) . This simple colour-based estima- 
tor over-predicts the mass-to-light ratio, particularly for younger, 
bluer galaxies. 



Figure [TU] shows the rest-frame K-h&nd mass-to-light 
ratio computed from the best-fit SEDs, as a function of the 
{g - r) colour for those WiggleZ galaxies with UKIDSS NIR 
photometry. Note, that the rest-frame Jf-band luminosity is 
calculated here using the best-fit SED for each galaxy and 
not the LEG template used to calculate the rest-frame FUV 
luminosity. The rest-frame Jf-band luminosity depends more 
critically on the exact choice of SED used to derive it. The 
iBell et all (|2003l ') colour-based estimator is also shown with 
an appropriate offset to match the IMF used in our SED fits. 
Although these mass-to-light ratios are computed using the 
KG04 code which does not include a prescription for dust 
emission in the NIR, we find in Table [T] that the KG04 de- 
rived stellar masses and MAGPHYS derived stellar masses 
are very similar. As MAGPHYS does include a prescription 
for dust emission in the NIR, we conclude that this emission 
has negligible impact on our calculated stellar masses and 
mass-to-light ratios. We find that even with the inclusion 
of secondary bursts, the stellar mass-to-light ratios derived 
from the SEDs, are lower than predicted by the colour-based 
estimator. The simple colour-based estimator over predicts 
the mass-to-light ratio by ~0.4 dex on average compared to 
the more sophisticated SED fitting approach in the case of 
these WiggleZ galaxies. Almost 75% of the WiggleZ galaxies 
have K-ha.nA mass-to- light ratios that are formally inconsis- 
tent with the simple colour-based estimator after taking into 
account the errors on M/Lk from the SED fitting. As seen 
in Figure fTOl this is particularly true for the younger, bluer 
galaxies. For the oldest and reddest WiggleZ galaxies with 
ages <;7Gyr, the SED derived M/Ljf and the colour-based 
estimator, agree reasonably well. As expected, the WiggleZ 
colour cuts select very few massive red galaxies that would 
populate the top right corner of Figure 1101 



5 MID-INFRARED LUMINOUS WIGGLEZ 
GALAXIES 

Before discussing the results of our SED fitting and stellar 
mass estimates within the broader context of galaxy evolu- 
tion, we briefly digress to look at the properties of the very 
small subset of WiggleZ galaxies that are found to be ex- 
tremely luminous at 12 and 22/im. There are 78 galaxies at 
0.3< z < 1.0 out of the total sample of 39,701, that are 
detected at >5a in WISE at 12 and 22^m. In Figure [TTl 
we plot their colours in the [ 3.4— 4.6] nm ve r sus [4 .6— 12]nm 
plane, taken from figure 12 of IWright eral] (|2010l '). We find 
unsurprisingly, that the mid infrared emission can be ac- 
counted for by the presence of an AGN in the majority of 
these galaxies with the WISE colours of the galaxies over- 
lappin g the QSO and Obscured AGN loci of IWright et al.l 
(|2010D . In order to select those galaxies where the AGN 
contamination should be less sig nificant, we apply a conser- 
vative cut of [3.4-4.6]pm < 0.7 (jStern et al.ll2012l l. which is 
also shown in Figure [TTl This leaves us with only 24 galaxies 
which may reasonably be assumed to be star-forming and 
lie in the starburst/LIRG regime of the colour-colour plane. 
We visually inspect the WISE images for these 24 galaxies 
and eliminate two which lie in the halos of bright sources in 
the WISE data. 

We fit the SEDs of the remaining 22 mid-IR bright 
star-formation dominated WiggleZ galaxies using the MAG- 
PHYS code which provides a consistent treatment of the 
UV, optical and IR emission in star-forming galaxies. Some 
example SED fits can be seen in Figure [T^] and encompass 
the range in properties seen in this small sample. In almost 
all cases, we find that the models under-predict the GALEX 
UV fiuxes of these galaxies and a UV excess is seen in the 
observed photometry relative to the models. This could be 
due to the effect of the Lya line on the GALEX fluxes of 
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[4.6-1 2]^_^ (Vega) 

Figure 11. The location of our WiggleZ mid-IR bright galax- 
ies in the [3.4— 4.6] mn versus [4.6— 12]^m colour-colour plane of 
jWright et al.|[201(J) . The WiggleZ galaxies are shown as the red 
circles. While most are found to be mid-IR bright on account of 
some AGN contamination, we isolate a population of 24 galaxies 
with [3.4-4.6] (.Lm < 0.7 that lie below the solid horizontal line in 
the figure and can reasonably be assumed to be star-formation 
dominated. 
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the galaxies or patchy dust extinction due to which some of 
the UV-hght still remains unobscured. 

The median stellar mass of these mid-IR bright WiggleZ 
galaxies is logio(M./Mo) = 10.8±0.6 inferred from SED fits 
from the UV through to the mid IR, and they typically 
have star formation rates in excess of lOOMgyr"^ and large 
dust extinction values. Their distribution in terms of redshift 
and observed UV-to-optical colours however, is the same as 
the rest of the WiggleZ sample. We conclude that a small 
fraction of the UV luminous WiggleZ galaxies at z ~0.7 are 
also infrared luminous with SEDs consistent with young, 
dusty starbursts where some of the young stellar population 
is still unobscured in the UV. This is consistent with what 
is found for the z ~ 1 Lyman Break Galaxy population 
some of which are also seen to be lumin ous in the Spitzer 
MIPS 24/im band (|Burgarella et al.ll2007l '). In Appendix B, 
we provide the positions, redshifts and WISE fluxes of these 
22 mid infra-red luminous WiggleZ galaxies along with the 
best-fit stellar masses and star formation rates derived using 
MAGPHYS. 
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6 DISCUSSION 

Having constrained the stellar masses of a large sample of 
UV-luminous galaxies at 0.3 < z < 1.0 and quantified the 
sensitivity of these mass estimates to assumptions made dur- 
ing the SED fitting process, we can now compare the stellar 
masses and SED fit parameters to other well-studied galaxy 
populations in order to place the WiggleZ galaxies within 
a global picture of galaxy evolution. Our study has concen- 
trated on galaxies that populate the most extreme end of 
the blue-cloud and, as such, represent a sample where simple 
SED models are not likely to provide a good representation 
of the galaxy physics. Despite this, we have shown the stel- 
lar masses to be extremely robust to different assumptions 



Figure 12. Best-fit SEDs from MAGPHYS for some of the Wig- 
gleZ galaxies that are very luminous at 12 and 22/xm. The plots 
show both the reddened and unreddened SEDs demonstrating 
that the galaxies have significant amounts of dust. Residuals from 
the reddened SED fits are plotted at the bottom of each panel. 



© 0000 RAS, MNRAS 000, 000-000 



Stellar Masses of WiggleZ Galaxies 19 



made during the SED fitting witli the median masses show- 
ing at most ~0.3 dex variation due to difTerences in the SPS 
models and /or inclusion of additional bursts of star forma- 
tion. Although the constraints on the star formation rates 
from the SED fitting are less robust with typical errors of 
>1 dex, the power of a large statistical sample such as ours 
is that the inferred median properties of the galaxies can 
reasonably be taken to provide a good representation of the 
sample as a whole. 

The median star formation rate derived using single 
component star formation histories and FAST is between 3- 
lOM0yr~^ regardless of the SPS model used and the number 
of photometric bands used in the SED fitting. The KG04 
code gives a median star formation rate of ~5 M0yr^^, 
which increases to only 5.3 M0yr~^ when secondary bursts 
are included. Finally, the subset of 1600 galaxies for which 
we compute best-fit SED parameters using MAGPHYS 
(Section I4.4.2p also has a median star formation rate of 6- 
7M0 yr~^ consistent with the results from the other SED 
fitting codes. 

In Figure [13] we plot the distribution of stellar masses 
and star formation rates for all ~40,000 WiggleZ galax- 
ies from this study. Spectroscopically confirmed Luminous 
Red Galaxies over the same redshift range represent a con- 
trasting population of extremely red galaxies for compari- 
son to the extremely blue WiggleZ galaxies. Stellar masses 
have been estima ted for large samp les of these selected 
from the 2SLAQ (iBanerii et a.l.ll201ol ') and SDSS-III-BOSS 



(|Maraston et al.l 20121 ') surveys. For comparison to the Wig- 
gleZ population, we choose LRGs from BOSS that are fit to 
have ongoing star formation (^40 per ce nt of the full BOSS 
LRG sample from iMaraston et al.ll201^ ) and also show the 
location of these in the Mt-SFR plane in Figure 1131 It is 
interesting to ask whether the very blue UV-luminous Wig- 
gleZ galaxies are the intermediate redshift analogues of the 
more distant Lyman Break and BX-selected galaxies at z~2 
fAdclbcrgcr ct al."2004!; IShaplev et ai]|2005l : lErb et al.ll2006l : 
[Ha bcrzc ttl ct al. 2011) as well as local UV-luminous galaxies 
(UVLGsV lHeckman et all (|2005h find that local UVLGs can 
be divided into large UVLGs and compact UVLGs. While 
the more massive large UVLGs have specific star formation 
rates (sSFR) sufficient to build their stellar mass over the 
Hubble time and therefore represent the most massive tail of 
star- forming disk galaxies like those in SDSS, the less mas- 
sive compact UVLGs have higher sSFRs and are typically 
observed in a starbursting phase. The local UVLGs have 
typical FUV luminosities of ~3xlO^°L0 while the higher 
redshift LBGs are slightly more luminous at FUV wave- 
lengths with typical FUV luminosities of '--^Gx 1O^''L0. Our 
sample of WiggleZ galaxies has a median FUV luminosity 
of ~3xlO^°L0 but at redshifts below 0.4, the median FUV 
luminosity is only ~9xlO^L0 whereas at z > 0.8 it is typ- 
ically ~7xlO^°L0. The low redshift WiggleZ galaxies are 
therefore less luminous than local UVLGs while the higher 
redshift WiggleZ galaxies have comparable FUV luminosi- 
ties to LBGs. How do the rest of their properties compare 
to the UVLG and LBG samples? 

In Figure 1131 we also plot the average ste l lar m asses 
and SFRs of UVLGs taken fromlHeckrnan et al.l (|2005|) and 
z~2 LBGs taken from Hcck man et al. l (|2005h : IShapiev et al.l 
(|2005h : lHabcrzcttl ct al.. (,20111 ') for comparison to the Wig- 
gleZ sample. These LBGs were selected either using GALEX 



(|Haberzettl et af] I2OIII) or the optic al colour-selection: 
BX/BM technique fShaplev et allbOOSh . The GALEX se- 
lected population is less biased against star-forming galax- 
ies with older stellar populations. Figure [131 also shows the 
evolving 'main-seq uence' of star-forming galaxie s taken from 
lElbaz et al.] (|200it ) at z = O. lNoeske etall ([20071 ') a.t z = 0.65 
roughly correspo nding to the median redshift of our Wig- 
gleZ sample, and iDaddi et al] ((2003) at 2 = 2. At higher 
redshifts, the main sequence systematically shifts to higher 
star formation rates for a given stellar mass, refiecting that 
high-redshift galaxies were more active than those in the 
local Universe. While it is difficult to empirically constrain 
any M,-SFR relation from our WiggleZ sample, which is 
highly incomplete in many regions of this parameter space, 
we can nevertheless ask where the WiggleZ galaxies lie rel- 
ative to the already derived M*-SFR relations from pre- 
vious works. We have already noted that large UVLGs lie 
on the main sequence at z — while the compact UVLGs 
show excess star formation relative to this main sequence. 
Similarly, the BX/BM LBGs also lie on the z = 2 main se- 
quence and therefore represent normal star-forming galaxies 
at these redshifts, while the GALEX selected LBGs show 
excess star formation relative to the z — 2 main sequence. 
Figure [13] shows that the WiggleZ galaxies represent a het- 
erogenous population with most lying at the upper end of 
the z = 0.65 main sequence, but also a cloud of galaxies 
above this main sequence, that are typically observed in a 
bursting phase (see also Jurek et al., in preparation). Al- 
though the scatter seen in Figure [TJ] is dominated by the 
scatter in the determination of the SFR from the SED fit- 
ting, there is evidence from those galaxies where the SFR is 
relatively well constrained, that WiggleZ galaxies with bluer 
optical (g — r) colours predominantly lie above the main se- 
quence. At a redshift of z = 0.7, a typical galaxy of 1O^''M0 
would need to be forming stars at a rate of ~1.4M0yr~^ in 
order to build up it's entire stellar mass through constant 
star formation over the age of the Universe. Most of the 
WiggleZ galaxies have star formation rates that are higher 
than this. 



We also note that the stellar masses of the Wig- 
gleZ galaxies are comparable to the local GALEX selected 
UVLGs, as well as the well-studied high-redshift LBG popu- 
lation. In Section [4.4.21 we found that the WiggleZ galaxies 
show a trend of increasing dust attenuation with increas- 
ing stellar mass and decreasing FUV luminosity as also ob- 
served in higher redshift UV-luminous galaxies. These Wig- 
gleZ galaxies may therefore reasonably be thought to repre- 
sent the intermediate redshift analogues of the local UVLG 
and LBG populations. 



We note by contrast that the star-forming LRGs at sim- 
ilar redshifts lie well away from the main-sequence and have 
star formation rates that are more than an order of magni- 
tude lower than the main-sequence of star-forming galaxies. 
This is consistent with these massive red galaxies having 
undergone the bulk of their stellar mass assembly at ear- 
lier epochs, when they presumably also experienced much 
higher levels of star formation. 
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Figure 13. Stellar mass versus star formation rate for all ^^40,000 WiggleZ g alaxies (left) as well as Luminous Red Galaxies (LRGs) 
in the BOSS survey with ongoing star formation fright: i Maraston et a lTM) . The greyscale re presents the den sity of points in this 
plane. These distributions are compare d to the averag e properti es of local UV-lum inous galaxies l|Heckman et al.ll2005[ ) as well as the 
2 ~ 2 Lyman break galajcy population llShaplev et al.|[2 005: Hab erzettl et al1l201lt). Dashed lines show the evolving main-sequence of 
star-forming galaxies at different redshifts jNoeske'eral ] i2007l: lElbaz et al.ll2007l : IPaddi et al.ll200Tt . The WiggleZ galaxies mostly lie at 
the upper end of the main-sequence at 2; ~ 0.7 with some clearly observed in a starburst phase. The LRGs on the other hand lie well 
below the main sequence. 



7 CONCLUSION 

We have conducted a detailed study of the stellar masses 
of a large sample of --^40,000 UV-Iuminous spectroscopically 
confirmed galaxies at0.3<z<1.0 selected within the WiggleZ 
survey. Around 30 per cent of the sample are matched to the 
wide-field NIR UKIDSS Large Area Survey and around 15 
per cent are additionally detected at 3.4 and 4.6^m in the 
all-sky WISE survey. The IR detected population represents 
the redder, more luminous and more massive end of the 
WiggleZ population with stellar masses that are on average 
0.6 dex larger for the UKIDSS detected galaxies and 0.8 
dex larger for the WISE detected galaxies. However, at the 
high redshift end of the WiggleZ sample at z > 0.7, there is 
evidence for a cloud of IR undetected galaxies that are just 
as massive and just as red in terms of their optical colours, 
as the IR detected galaxies. In addition, we find a small 
sample of 22 galaxies which are also extremely luminous 
at 12 and 22^m with SEDs consistent with dusty starburst 
galaxies where some of the younger stellar population is still 
unobscured in the UV. 

As the WiggleZ galaxies represent the most extreme 
end of the blue cloud population at these redshifts, where 
SED fitting is likely to be the most problematic, we quantify 
the sensitivity of our stellar mass estimates to assumptions 
made during the SED fitting process. In particular we find 
that: 

• The effect that the NIR photometry has in constrain- 
ing the stellar masses, depends on the priors assumed in 
the SED fitting. With the SED fitting code, FAST used in 
conjunction with the BC03 and Maraston models, we find 
that the stellar mass constraints are improved on addition 
of the NIR photometry as this allows us to tighten the up- 
per bound on the stellar mass estimates. With the KG04 
code and the PEGASE.2 SPS models, we find that the NIR 
data makes little difference to the stellar masses and their 
corresponding errors. Regardless of the choice of SED fit- 



ting code, the optical and optical-|-NIR derived masses are 
consistent within the la errors for >68% of the galaxies as 
expected. With the inclusion of NIR photometry, the mass 
estimates from both FAST and KG04 agree very well. 

• The addition of UV photometry from GALEX makes 
very little difference to both the median stellar masses and 
the median la errors on these stellar masses, even for these 
extremely blue WiggleZ galaxies selected in the UV. How- 
ever, there is a small population of IR-detected WiggleZ 
galaxies where the lack of UV photometry leads to best-fit 
SEDs with considerably higher star formation rates, dust ex- 
tinctions and younger ages than when the UV photometry 
is included in the fitting. A consequence of fitting younger 
stellar populations without the UV data, is that the stellar 
masses of these galaxies are also under-estimated without 
the UV. 

• The choice of SPS model can affect the stellar mass 
estimates by ~0.3 dex. The Maraston models which include 
TP-AGB stars result in stellar masses that are ~0. 1-0.2 dex 
lower than those inferred using the BC03 models with the 
differences largest for the more FUV luminous subset of 
WiggleZ galaxies, which are typically older. The PEGASE.2 
models produce stellar masses that are 0.2(0.3) dex higher 
than BC03 (Maraston). The effect is more pronounced for 
galaxies with fainter FUV luminosities. 

• The inclusion of nebular emission in the SED models, 
lowers the stellar masses by only ^^0.02 dex on average. This 
offset is smaller than that seen in the case of high redshift 
LBGs. 

• The choice of IMF can affect the stellar mass by up 
to 0.3 dex. The Salpeter IMF produces the largest stellar 
masses. The Kroupa IMF results in stellar masses that are 
^0.1-0.2 dex lower than Salpeter depending on the SPS 
model being fit. The BG03 and Chabrier IMFs produce very 
similar stellar masses with a median difference of <0.001 
dex between them. These masses are ~0.24 dex lower than 
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those derived using the Salpeter IMF and ~0.1 dex lower 
than those derived using the Kroupa IMF. 

• For single component SFHs, the stellar masses are in- 
sensitive to the choice of star formation history although 
we find as expected that the truncated SFH does not pro- 
vide good values for these actively star-forming galax- 
ies. The addition of bursts of star formation on top of the 
smooth underlying SFH however, results in stellar masses 
that are up to 0.3 dex larger. The effect is more pronounced 
at higher FUV luminosities where secondary bursts can ef- 
fectively hide an older, more massive and more evolved stel- 
lar population. When bursts are allowed, the typical burst 
fraction is estimated to be ~1 per cent of the total stellar 
mass in these galaxies. 

• Although the dust extinction parameter. Ay is not as 
well constrained as the stellar mass, we note a trend of in- 
creasing dust attenuation with increasing stellar mass and 
decreasing FUV luminosity in these galaxies. These trends 
are consistent with what has been found in higher redshift 
UV-luminous galaxies. 

• We compare our best-fit rest-frame K-h&nA mass-to- 
light ratios from the SED fits to the predictions from 
simple op t ical c olour based estimators such as those of 
iBell et all (|2003h . We find that the colour-based estimator 
over-predicts the mass-to-light ratio for ~75% of the Wig- 
gleZ galaxies where the results from this estimator and the 
SED fits are formally inconsistent given the Icr errors. The 
median difference in M/L/f between the simple colour-based 
estimator, and the more sophisticated SED fitting approach, 
is ~0.4 dex. The inconsistencies are most pronounced for the 
bluer WiggleZ galaxies with young best-fit ages and galaxies 
with ages >7Gyr generally have M/Lk consistent with the 
colour-based estimator. 



We conclude that our stellar mass estimates have typ- 
ical dispersions of ~0.1-0.3 dex as a result of changes in 
the input parameters of the SED fitting and they are there- 
fore extremely robust to these changes even for these ex- 
tremely blue galaxies. While the star formation rates have 
much poorer constraints from the photometric data, the ad- 
vantage of a large sample such as ours is that the median 
values inferred from the SED fits, should provide a reason- 
able representation of the average properties of the WiggleZ 
sample as a whole. We find that the WiggleZ galaxies have 
star formation rates of ~3-lOM0yr~^ and these estimates 
are consistent between the different SED fitting codes tested 
in this work. Although these star formation rate estimates 
have very large errors for individual galaxies, we compare the 
distribution of the WiggleZ galaxies in the M*-SFR plane to 
other well-studied samples such as Luminous Red Galaxies 
over the same redshift range, local UV-luminous galaxies, 
and more distant Lyman Break Galaxies. We find that the 
WiggleZ population, on average lies at the upper end of the 
main-sequence of star- forming galaxies at z ~ 0.7. However, 
there is evidence for some of the bluer galaxies lying well 
above the main sequence and therefore being observed in 
a starburst phase (Jurek et al., in preparation). The stel- 
lar masses of our WiggleZ galaxies are comparable to both 
2 ~ 2 LBGs and compact UV-luminous galaxies in the lo- 
cal Universe suggesting that this population can be taken to 
represent a reasonable intermediate redshift analogue that 



straddles the redshift space between the well studied local 
UVLG and distant LBG populations. 

We conclude that the combination of current cosmolog- 
ical volume spectroscopic surveys and wide-field photomet- 
ric surveys, provides a promising route for constraining the 
physical properties of very large, statistically robust samples 
of intermediate redshift galaxies. Understanding the prop- 
erties of these z~l galaxies is essential for bridging the gap 
between the local and high redshift Universe, as well as de- 
vising effective target selection criteria for the next genera- 
tion of wide-field spectroscopic surveys. 

The stellar masses in this paper will be made publicly 
available with the next WiggleZ data release. 
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APPENDIX A: EXAMPLE SED FITS 

For illustrative purposes, we present some examples of the 
SED fits for the WiggleZ galaxies derived using the KG04 
code with multi-component SFHs and the PEGASE.2 SPS 
models including nebular emission. These example SEDs 
have specifically been chosen to encompass the full range 
in properties of our sample in terms of redshift, colour, stel- 
lar mass and age and also include some examples of poorly 
fit SEDs. 



APPENDIX B: MID-INFRARED LUMINOUS 
WIGGLEZ GALAXIES 

We presented a small sample of mid infrared luminous Wig- 
gleZ galaxies detected at 12 and 22/^m in WISE. Based on 
the WISE colours, these galaxies do not seem to have a 
significant AGN component and can be considered to be 
star-formation dominated. Here we provide a table of the 
positions, redshifts and WISE fluxes for these 22 galaxies in 
Table Bl. The best-fit stellar mass and star formation rate 
derived using MAGPHYS are also given in the same Table. 
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Figure Al. Example SED fits encompassing tiie full range of properties in our sample and showing some examples of galaxies that are 
poorly fit as well. These fits are derived using the KG04 code and the PEGASE.2 SPS models and including nebular emission lines. 
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Table Bl. Summary of WISE fluxes, best-flt stellar mass and star formation rates 
22/im in WISE, and where the colours rule out a significant AGN component. 
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